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We demonstrate that first-order phase transitions in 1+1-dimensional nonequilibrium systems with 
fluctuating ordered phases are impossible, provided that there are no additional conservation laws, 
long-range interactions, macroscopic currents, or special boundary conditions. Since minority islands 
in the ordered phase of such systems can only shrink by short-range interactions, it is impossible to 
stabilize a fluctuating ordered state. The apparent first-order behavior turns out to be a transient 
phenomenon, crossing over to a continuous transition after very long time. As examples we consider 
the triplet creation model 3X — » AX, X — » 0, the annihilation/fission process 2X — > 3X, 2X — > 0, 
as well as spreading on a diffusing background X+Y — > 2Y, Y — > X. 



I. INTRODUCTION 

In equilibrium statistical mechanics it is well known 
that systems undergoing a first-order phase transition in 
high space dimensions d may exhibit a second-order tran- 
sition below their upper critical dimension. The reason 
is that in low dimensions fluctuations become more im- 
portant and may destabilize the ordered phase. A very 
similar situation emerges in the case of nonequilibrium 
phase transitions 0]]. In this context the question arises 
under which conditions first-order phase transitions can 
be observed in one spatial dimension. The purpose of 
the paper is to point out that various 1+1-dimensional 
nonequilibrium models, which were believed to exhibit a 
first-order transition, cross over to a continuous transi- 
tion after very long time. 

Dynamic random processes with first-order phase tran- 
sitions are usually characterized by several stable or- 
dered phases. For example, in the subcritical regime 
< T < T c of the two-dimensional kinetic Ising model 
there are two stable magnetized states. To ensure their 
stability, the Ising model provides a robust mechanism 
eliminating islands of the minority phase generated by 
thermal fluctuations. This mechanism relies on the fact 
that the boundary of an island costs energy, leading to 
an effective surface tension. Attempting to minimize its 
energy, the island is subjected to an attracting 'force' 
and begins to shrink. It is important to note that this 
long-range force decreases algebraically as 1/r with the 
typical radius r of the island so that thermal fluctua- 
tions of any size are safely eliminated. Because of the 
^-symmetry under spin reversal, both ordered phases 
are equally attracting. Thus, starting from a disordered 
state with zero magnetization, we observe coarsening pat- 
terns of ordered domains. However, if an external field 
is applied, one type of minority islands becomes unsta- 
ble above a certain critical size. Since there is a finite 
probability to generate such islands by fluctuations, one 
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FIG. 1. Schematic phase diagram of Toom's 

north-east-center voting model. 

of the two ordered phases eventually takes over, i.e., the 
system undergoes a first-order phase transition. 

Turning to nonequilibrium systems, the mechanism for 
the elimination of minority islands may be even more ro- 
bust. An interesting example is Toom's two-dimensional 
north-east-center voting model ||||]. Again there are 
two stable ordered phases. In contrast to the kinetic 
Ising model, where straight domain walls perform an un- 
biased diffusive motion, the interfaces in Toom's model 
propagate in a preferred direction, wherefore the process 
is out of equilibrium. As a key property of the model, 
this propagation velocity depends on the orientation of 
the domain walls. Using this anisotropy, the dynamic 
rules of the model are designed in a way that minority 
islands quickly assume a triangular shape and begin to 
shrink with constant velocity. Thus, the effective 'force' 
by which an island shrinks is independent of r, leading to 
much more stable phases as in the standard kinetic Ising 
model. Consequently, the ordered phases remain stable 
even if an oppositely oriented external field H is applied. 
In order to flip the whole system, the intensity of the 
field has to exceed a certain critical threshold. The cor- 
responding phase stability boundaries (also called spin- 
odal lines) are sketched in Fig. [p. Between these lines, 
the two ordered states are both thermodynamically sta- 
ble, i.e., they coexist in a whole region of the parameter 
space. Crossing the coexistence region by varying the 



1 



external field, the magnetization of the system follows in 
a hysteresis loop. Similar phenomena can be observed in 
certain models for nonequilibrium wetting |I| . 

Obviously, both mechanisms for the elimination of mi- 
nority islands - surface tension and anisotropic propa- 
gation velocities - can only be implemented in at least 
two spatial dimensions. Therefore, it is interesting to 
investigate the question under which conditions first- 
order phase transitions can be observed in one spatial di- 
mension. For example, the 1+1-dimensional Ziff-Gulari- 
Barshad model || for heterogeneous catalysis is known to 
exhibit a first-order phase transition which relies on the 
interplay of three different kinds of particles. Similarly, 
a recently introduced model for phase separation on a 
ring jq] uses three different species of particles. Another 
example is the so-called bridge model [[?J for bidirectional 
traffic on a single lane, where special boundary conditions 
induce a discontinuous transition in the currents. Even 
more subtle is the mechanism in the two-species model 
introduced in a recent paper by Oerding et al, where a 
first-order phase transition is induced by fluctuations || . 
Therefore, in attempting to comprehend the full range 
of first-order phase transition under nonequilibrium con- 
ditions, it would be interesting to seek for the simplest 
1+1-dimensional model which exhibits a discontinuous 
transition. By 'simple' we mean that such a model should 
involve only one species of particles evolving by local 
dynamic rules without macroscopic currents, conserva- 
tion laws, and unconventional symmetries. Moreover, the 
choice of the boundary conditions should be irrelevant. 

The prototype of such a dynamic process is the one- 
dimensional Glauber-Ising model at zero temperature in 
a magnetic field. This model, which is also referred to 
as compact directed percolation (CDP), can be defined 
as follows. Sites of a one-dimensional lattice can be in 
two different states T and J,. For each update a pair of 
adjacent sites is randomly selected. If the two spins are 
in opposite states, the are aligned with the probabilities 
p and 1 — p: 

TUT TT, TUT ^ II • (l) 

Obviously, the parameter H = p—l/2 plays the role of an 
external field. Since there is no temperature, the model 
has two absorbing configurations, namely the fully mag- 
netized states A = ... |J,| ... and B = ... TTT ••■ ■ Once 
the system has reached one of the two absorbing config- 
urations, it is trapped and will remain there forever. 

In the off-critical case H ^ 0, one of the two absorbing 
states is stable while the other one is unstable against 
small perturbations. For example, for H < 0, S-islands 
tend to grow while A-islands shrink. Thus, starting with 
random initial conditions, the system approaches the 
thermodynamically stable state B in an exponentially 
short time. In the vicinity of the phase transition this 
time scale diverges as \H\~ . Right at the critical point, 
the two absorbing states of the Glauber model arc only 



marginally stable against perturbations. For example, by 
flipping a single spin in a fully magnetized domain, we 
create a pair of kinks. These kinks perform an unbiased 
random walk until they annihilate one another. Thus, 
minority islands do not shrink by virtue of an attract- 
ing force, rather they are eliminated solely because of 
the fact that random walkers in one dimension always 
return to their origin. Consequently, the lifetimes r of 
minority islands are finite and distributed algebraically 
as P(t) ~ t -1 / 2 . It is important to note that this mech- 
anism allows islands to reach a macroscopic size of the 
order y/r during their temporal evolution. Therefore, av- 
eraging over many independent samples, the mean size 
of such minority islands approaches a constant value. As 
a consequence, the ordered phase is only marginally sta- 
ble against perturbations. In fact, introducing a small 
rate for spontaneous spin flips, the first-order transition 
in the Glauber-Ising model is lost. The same happens in 
certain nonequilibrium Ising models with two absorbing 
states subjected to an external field ||,[l0) ■ 

Is it possible to observe first-order transitions in fluc- 
tuating 1+1-dimensional two-state systems? Obviously, 
such a model requires a much more robust mechanism 
for the elimination of minority islands. A simple ran- 
dom walk of a pair of kinks is not sufficient, rather there 
has to be an attracting force which prevents small island 
from growing. It seems that such a mechanism is difficult 
to implement. In fact, various models, which have been 
suggested in the past, display a discontinuous transition 
only in d > 2 dimensions jl|^,|||] . 

A frequently cited exception is the 1+1-dimensional 
triplet creation process introduced by Dickman and 
Tome @: 

3A -> AX , X -> . (2) 

In this model the high-density phase is not strictly ab- 
sorbing, rather islands of unoccupied sites are sponta- 
neously created in the bulk so that the active state fluc- 
tuates. In numerical simulations it was observed that 
above a certain tricritical point the second-order phase 
transition line splits up into two spinodal lines, where the 
transition becomes first order. Moreover, the order pa- 
rameter seemed to follow a hysteresis loop when the pa- 
rameter for offspring production was varied. Apparently, 
the highly nonlinear particle creation process 3X — > AX 
is able to stabilize the high-density phase, eliminating is- 
lands of the minority phase. Similarly, Carlon et al. 
reported a first order transition in the 1+1-dimensional 
annihilation/fission process on a lattice fl8|| 

2X -> 3A , 2X -> . (3) 

For low values of the diffusion constant they found a con- 
tinuous transitions which becomes first order for high dif- 
fusion rates above a certain tricritical point. 

In the present work we argue that this type of first- 
order behavior in 1+1-dimcnsional models is only a tran- 
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sient phenomenon. Although mean field approximations 
predict a discontinuous behavior, the transition crosses 
over to a continuous transition for any value of the pa- 
rameters (i.e., for any value of the diffusion constant in 
the examples mentioned before). 

The paper is organized as follows. Discussing the dy- 
namics of interacting domain walls in Sect, uA we ar- 
gue that short-range forces are not sufficient to stabilize 
fluctuating ordered phases in one spatial dimension. In 
Sees. [II and [V we revisit the triplet creation process and 
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the annihilation/fission process, both being examples of 
models with short-range forces. As a third example a 
spreading process on a diffusing background will be dis- 
cussed in Sec. [v|. Although in this case long-range forces 
lead to a coarsening process, the transition turns out to 
be continuous as well. 



II. INTERACTING DOMAIN WALLS IN ONE 
DIMENSION 



Short-range forces between domain walls 

To explain the impossibility of first-order transitions in 
the processes (|])-(||) from a phenomenological point of 
view, let us consider a hypothetical one-dimensional sys- 
tem with two ordered phases A and B. Without loss of 
generality we assume the A-phase to be absorbing while 
the B-phase fluctuates, i.e. small A-islands are spon- 
taneously created in the bulk of the B-phasc. Further- 
more, we assume that there is a robust mechanism which 
eliminates minority islands generated by fluctuations, en- 
suring the stability of phase B. For a system with these 
properties, let us consider an initial state where half of 
the system is in phase A while the other half is in phase B 
(see Fig. |J). Both phases are separated by a domain wall. 
In contrast to the Glauber-Ising model, this domain wall 
is not necessarily associated with a single broken bond, 
rather is may be 'smeared out' over a certain region in 
space. However, since both phases are attracting, the 
domain wall remains localized, i.e., it has a finite width. 
Consequently, the derivative of the order parameter pro- 
file has exponentially decreasing tails. Obviously, the 
critical point of such a model corresponds to a situation 
where the domain wall diffuses without bias. 

In one dimension a minority islands of phase A inside 
phase B can be considered as a pair of such domain walls. 
In the Glauber-Ising model the domain walls do not in- 
teract unless they annihilate at the same site. In our 
hypothetical model, however, the domain walls do inter- 
act, giving rise to an attracting 'force' eliminating islands 
of the minority phase. However, since the domain walls 
have a finite width, we expect this 'force' to have a finite 
range ro . More precisely, for large island sizes r we expect 
the strength of the interaction to decrease exponentially 
as exp(— r/rp) or even faster. 
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FIG. 2. Schematic profile of the coarse-grained particle 
density p(x) and its derivative p'(x) at a domain wall between 
phase A and phase B in a one-dimensional system. 



Instability of fluctuating ordered domains 

The essential problem arises precisely at this point. In 
one spatial dimension an exponentially decreasing short- 
range force between two domain walls is not sufficient 
to stabilize a fluctuating ordered phase. In order to un- 
derstand this point, let us consider an island of size r 
which grows by one step with rate 1 — aexp(— r/fo) and 
shrinks with rate l + aexp(— r/r^). In this toy model the 
parameters a and tq control the strength and the range 
of the force, respectively. The time evolution of the cor- 
responding probability distribution P r (t) is given by the 
equation)^ 

^P r (t) = -2P r (t) + (\ + ae-^l^)P r+1 {t) 
+ (l-ae- ( - r - 1 V ro )P r _ 1 (t) 

with the boundary condition Po(t) ~ 0. Starting with 
a single spin flip P r (0) — <5 r ,i, we obtain a broadening 
distribution of island sizes. Without a short range force 
(i.e. ro = 0), Eq. (^) reduces to the discrete diffusion 
equation with a Dirichlet boundary condition at r = 0. 
In this case the asymptotic solution reads 



Pr(t) 



Ct- 3/2 cxp(-r 2 /4i) : 



(5) 



where C is a certain amplitude factor. In presence of a 
short-range force with ro > 0, this asymptotic solution 
remains valid - the only change is in the prefactor C. 
Although C decreases exponentially with increasing ro, 
it is always positive. Therefore, there is a small but finite 
probability that the two domain walls become effectively 
unbound, producing a macroscopic minority island with 
a size of the order y/r. Thus, the ordered phase is only 
marginally stable, irrespective of the interaction range ro. 
As a consequence, the B-phase eventually disintegrates, 
approaching the absorbing state A. This contradiction in 
our hypothetical model demonstrates that in one spatial 
dimension it is impossible to stabilize fluctuating ordered 



*In principle there may be two different length scales on both 
sides of the domain wall. 
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phases by short-range interactions. It should be empha- 
sized that these arguments are not valid in higher dimen- 
sions where the dynamics of domain walls may depend 
on the local curvature. 

Properties of the continuous transition 

Since absorbing islands may reach a macroscopic size, 
our hypothetical model approaches the absorbing state. 
In order to prevent the system from reaching the absorb- 
ing state, we may now slightly modify the parameters 
such that B domains tend to grow. Since the A-phase 
is absorbing, the resulting competition between growing 
B-domains and spontaneously created A-islands suggest 
a crossover to directed percolation (DP) |p^| , po| , which is 
the generic universality class of continuous phase tran- 
sition into absorbing states. For instance, as will be 
shown in Sec. Ill, the triplet creation process displays 



such a crossover to DP. In general, the universality class 
of the continuous transition will depend on the symme- 
try properties and the conservation laws of the process 
under consideration. As an example, we will discuss the 
annihilation/fission (see Sec. IV) which crosses over to a 
non-DP transition. 

Numerical checks 

Obviously these considerations are only valid if the or- 
dered states are sufficiently attracting to confine domain 
walls to a finite region with a typical size r$. There are 
several possibilities to verify this assumption. The sim- 
plest method would be to perform high-precision Monte 
Carlo simulations in order to confirm the crossover to 
a continuous transition. Alternatively, we may check 
whether a compact £?-phase disintegrates as predicted. 
As a more sensitive check, one could also measure the 
distribution of the sizes of minority islands. Since spin 
flips occur spontaneously in the ordered phase, we expect 
this probability distribution to become stationary in the 
limit t — > oo. Solving the stationary problem of Eq. (Q), 
we obtain the distribution 



r-1 

n.ir 



-r'/ro 



-(r' + l)/r 



Pi 



(6) 



This distribution decays exponentially for small r until 
it saturates at a small but finite constant in the limit 
r — ► co, expressing the fact that minority islands of all 
sizes will be formed. A distribution of this form observed 
in a numerical simulation can be considered as a hallmark 
for the presence of a short-range force between domain 
walls. 



III. THE TRIPLET CREATION PROCESS 

3X -> AX, X -> 

We are now going to apply these tests to the triplet cre- 
ation model introduced by Dickman and Tome |lq ]. It 
evolves by random sequential updates and is controlled 
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FIG. 3. Numerical simulation of the triplet creation pro- 
cess. Starting with an island of 20 particles the survival prob- 
ability P(t) (left) and the mean number of particles N(t) 
(right) of the cluster are averaged over 300 independent runs. 
The data points are multiplied with powers of t in a way that 
DP critical behavior corresponds to horizontal lines (see text). 



by two parameters, namely a rate for offspring produc- 
tion A and a diffusion constant D. For D < 0.85 the 
transition was found to belong to the universality class 
of directed percolation, whereas for D > 0.85 a first order 
transition was reported. Moreover, the phase transition 
line seemed to split up into two spinodal lines X±(D). In 
the following we show that the transition crosses over to 
DP after very long time. We restrict the analysis to the 
case D = 0.9, where Dickman and Tome reported the 
stability limits A_ = 10.12(1) and A+ = 10.30(2). 

As a first test, we perform Monte Carlo simulations 
starting with a single island of particles at the origin [^l] . 
It turns out that it is necessary to start with an island 
of several particles since otherwise the survival probabil- 
ity of the cluster is extremely small. At criticality the 
survival probability P(t) and the average number of par- 
ticles N(t) are expected to obey asymptotic power laws 
of the form 



Pit) ~ P 



N{t) ~ i" 



(7) 



with 8 — 1/2, r\ = for CDP (indicating a first-order 
transition) and 6 ~ 0.159, rj ~ 0.313 in the case of a DP 
transition, respectively. The results of the simulations 
are shown in Fig. ||[ Obviously, there are three different 
temporal regimes. In the first 100 time steps, the island 
survives with certainty due to the large initial size of 20 
sites, followed by a rapid decrease of P(t) faster than 
X/yt. This part of the temporal evolution is expected 
to be non-universal. Then the model enters a second 
regime extending over two decades from 10 3 to 10 5 time 
steps, where P{t) decays as 1/yt and N(t) stays almost 
constant. This is the time window where the model be- 
haves essentially in the same way as a zero-temperature 
Glauber-Ising model so that the transition appears to 
be discontinuous. In the following two decades from 10 5 
to 10 7 time steps, however, we observe a slow crossover 
to DP exponents. Our estimate for the critical point 
A = 10.145(10) lies between the stability limits A± re- 
ported in [Q. With the computer technology ten years 
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ago, Dickman and Tome could only go up to 22000 time 
steps so that they were unable to observe the crossover to 
DP. A similar crossover phenomenon from an initial tran- 
sient over CDP to DP was recently observed in certain 
models for flowing sand on an inclined plane |22] . 

As a second check, we demonstrate that the fluctuating 
phase disintegrates, generating the typical patterns of DP 
clusters after very long time. To this end we introduce 
a novel type of space-time plots which can be used to 
visualize the scaling properties of critical clusters in sys- 
tems with absorbing states (see Fig. ||). Starting with a 
localized seed (an island of 20 particles) at the origin and 
simulating the process up to 10 7 time steps, the rescaled 
position of the particles x/t 1 ^ is plotted against log 10 t, 
where z = v\\ /v± is the dynamic exponent of the pro- 
cess under consideration. Here we choose z = 2 since the 
domain walls are expected to diffuse. By rescaling the 
spatial coordinate x, the cluster evolves within a strip 
of finite width. Unlike linear space-time plots the scale- 
invariant representation of a cluster allows one to survey 
more than six decades of the temporal evolution. As can 
be seen, the figure is consistent with the Monte Carlo 
simulation results. After an initial transient the cluster 
appears to be compact, behaving essentially in the same 
way as in a Glauber-Ising model at zero temperature. 
The boundaries, which are smeared out by diffusion, be- 
come sharper as time proceeds. The first minority island 
with a macroscopic size appears after 2 • 10 4 time steps, 
where the crossover to DP begins. However, the typical 
patterns of a DP cluster can only be seen after 10 6 time 
steps. 

As a third check, we measure the stationary distribu- 
tion P r of island sizes in order to verify the predictions 
of the toy model in Sect. ||. Starting with a fully oc- 
cupied lattice, the sizes of inactive islands are measured 
in a critical system of 10000 sites after an 'equilibration' 
time of 10 6 time steps. The result is shown in Fig. |^ as a 
solid line. As can be seen, the distribution P r decays ex- 
ponentially or even faster by almost five decades until it 
crosses over to a slowly decreasing function. Fitting the 
parameters a and tq, this curve may be compared with 
the stationary solution of the toy model in Eq. (^]). Al- 
though the curves do not coincide, they show the same 
qualitative behavior. This supports the point of view 
that domain walls in the 1+1-dimensional triplet creation 
process attract one another by a short-range force so that 
the fluctuating active phase cannot be stabilized. 



IV. THE ANNIHILATION/FISSION PROCESS 

2X -> 3X, 2X -> 

Some time ago Howard and Tauber introduced the so- 
called annihilation/fission process which exhibits an un- 
conventional nonequilibrium phase transition . Their 
original motivation was to consider a reaction-diffusion 
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FIG. 4. The triplet creation process 3X -> 4X, X -> 
at the critical point. The figure shows the temporal evolu- 
tion over seven decades starting with a compact island of 20 
particles in the center. The graphs shows the rescaled posi- 
tion x/\/t of the particles as a function of log(t) measured in 
a single run. Up to 10 4 time steps the cluster has the form 
of a compact diffusing cloud of particles with very small is- 
lands of unoccupied sites generated by fluctuations. The first 
macroscopic minority island emerges after 2 • 10 4 time steps, 
indicating the beginning of the crossover to DP which extends 
up to 10 6 steps. Only in the last decade, where the thickness 
of active branches is small compared to the lateral cluster 
size, we observe the typical patterns of a directed percolation 
process. 
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FIG. 5. Distribution of the sizes of inactive islands in the 
critical triplet creation model (solid line), compared to a fit 
of Eq. (^) (dashed line). 



process which allows one to interpolate between 'real' and 
'imaginary' noise in the corresponding Langevin equa- 
tion. Performing a field-theoretic renormalization group 
analysis they predicted non-DP critical behavior at the 
transition. 

More recently, Carlon et al. |L7| studied a lattice ver- 
sion of the model which is defined by the dynamic rules 

AA0, 0AA — > AAA with rate (1 - p) (1 - d)/2 

AA -> 00 with ratep(l - d) (8) 
A0 <-> A with rate D . 

Performing a density matrix renormalization group anal- 
ysis they arrived at the conclusions that for low values 
of the diffusion constant D the transition is continuous 
while it becomes first order for higher values of D above 
a certain tricritical point. The paper by Carlon et al. 
released a debate concerning the universality class of the 
continuous phase transition. Because of a numerical co- 
incidence of two out of four critical exponents, the au- 
thors concluded that the transition should belong to the 
parity-conserving universality class ]23|-|28|]. This result 
was questioned in Ref. p9| since there is neither a parity 
conservation law nor a ^-symmetry in the AF process. 
Moreover it was pointed out that a first order transi- 
tion might not exist in one spatial dimension. Increasing 
the diffusion constant, continuously varying critical ex- 
ponents were reported in |5(j| . For large values of D these 
exponents seemed to approach certain values which can 
also observed in cyclically coupled DP and annihilation 
processes |3l]]. So far there is no numerical evidence for 
a discontinuous transition. 

It should be emphasized that the AF process has two 
non-symmetric absorbing states, namely the empty lat- 
tice and the state with a single diffusing particle. The 
unusual critical behavior is related to the fact that soli- 
tary particles may diffuse over large distances before they 
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FIG. 6. Scale-invariant plot of a space-time history of the 
annihilation/fission process for D = 0.9 at the critical point 
p c = 0.233(2). As can be seen, the high density phase disin- 
tegrates, creating macroscopic islands of the minority phase 
in the bulk. 

meet another particle in order to annihilate or to release 
a new avalanche. Therefore, the absorbing state may be 
'less attracting' than in the triplet creation model. Be- 
cause of these special properties, it is not immediately 
clear whether the assumptions of Sect, [n] apply. For ex- 
ample, solitary diffusing particles could lead to an effec- 
tive long-range interaction between domain walls. How- 
ever, as demonstrated in Fig. |^, even for high values of the 
diffusion rate a compact island quickly disintegrates, gen- 
erating nontrivial patterns of active patches and solitary 
particles. This supports the claim by Odor |3(| that there 
is no first-order phase transition in the AF process for any 
< D < 1. The question whether the observed contin- 
uously varying critical exponents are related to asymp- 
totically well-defined properties or simply to crossover ef- 
fects between an an effective first-order behavior on short 
scales and an asymptotic continuous transition remains 
open. 
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V. SPREADING PROCESS ON A DIFFUSING 
BACKGROUND 



Let us finally investigate the following two-species par- 
ticle process with random-sequential updates: 

X + Y -c Y + Y with rate a 

Y -> X with rate j3 (9) 

X + <-> + X with rate L» x 

y + <-> + y with rate D y 

This reaction diffusion scheme may be interpreted as 
a spreading process of Y particles on a diffusing back- 
ground of X particles. In this process the total number of 
particles is conserved. Therefore, the spreading proper- 
ties are controlled by the density of particles in the initial 
state p = px + py, while the density py plays the role of 
an order parameter. For low values of p the particles are 
sparsely distributed so that the system quickly evolves 
into an 'absorbing' state without Y particles. For high 
values of p the spreading process X + Y — > 2Y dominates, 
leading to a stationary density py > in a sufficiently 
large system. Both phases are separated by a nonequi- 
librium phase transition at a certain critical density p c 
which depends on the parameters a, (3 and the diffusion 
rates Dx and Dy. 

The critical behavior of the model depends signifi- 
cantly on the ratio of the diffusion rates Dx and Dy. 
If both rates are equal, the particles diffuse in the same 
way as in an ordinary exclusion process, whereas the la- 
bels X and y can be considered as 'colors' with separate 
dynamic rules. This special case has been considered in 
Ref. , where a continuous phase transition was found. 
Recently, van Wijland et al. were able to confirm these 
results by a field-theoretic analysis |3^]. Moreover, they 
investigated the general case Dx ^ Dy, For Dx < Dy 
they predicted a continuous phase transition with non- 
DP critical exponents, while for Dx > Dy it turned out 
to be impossible to perform a controlled field-theoretic 
analysis. For this case van Wijland et al. conjectured 
that the transition might be discontinuous. 

Obviously, the arguments used in Sect, [n] cannot be 
applied to this model since the particle conservation law 
introduces effective long-range interactions. Thus, the 
model may well be a candidate for a first-order transition 
even in one spatial dimension. To investigate this ques- 
tion in more detail, we analyzed the two-species model 
numerically. Choosing a = 0.2, (3 = 1 and Dy = 1 — Dx, 
the process (||) is controlled by two parameters Dx — D 
and p. The Y particles may be considered as the active 
sites of a spreading process running on top of a diffusing 
background of X particles. For D < 0.5 the transition 
is characterized by fractal clusters reminding of DP, but 
apparently in a different universality class. For D > 0.5, 
however, there is indeed a clear signature of first-order 




FIG. 7. Semi-logarithmic plot of the spreading process (||) 
with 1000 sites for D = 0.9 in the active phase p — p c . The 
figure shows the dynamics of the Y particles, while the diffus- 
ing background of X particles is not visible. As can be seen, 
islands of Y particles coarsen until a single island is left. 



behavior. As shown in the space-time plot in Fig. [j], the 
y particles are confined to certain vertical 'channels'. As 
time proceeds, small channels disappear while larger ones 
tend to grow, leading to a coarsening process. 

The mechanism behind this coarsening process can be 
understood as follows. Because of Dy < Dx, a high 
concentration of Y particles reduces the effective diffu- 
sion rate. Since a locally reduced diffusion rate leads to 
particle clogging, the local density of particles increases. 
Thus there is positive feedback amplifying the particle 
density (and hence the effective spreading rate) in regions 
where the spreading process is active. Because of this in- 
stability, the y particles tend to form compact domains 
with a high particle density. In d > 2 spatial dimensions, 
such an instability may lead to a first-order transition. 
In one-dimensional systems, however, the nonlinear am- 
plification mechanism alone is not sufficient, especially 
when the high-density phase fluctuates. Obviously, the 
essential mechanism behind the coarsening process is the 
conservation of the total number of particles. For exam- 
ple, let us assume that a high-density domain has been 
created by fluctuations. Because of the instability, the 
spreading process is supercritical inside the domain. Par- 
ticle conservation ensures that such an island cannot split 
up into two separate islands since there is no way to re- 
duce the density inside. 
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FIG. 8. Spreading process on a diffusing background. 
Stationary density of Y-particles as a function of the to- 
tal particle density p for various values of the diffusion rate 
D — D x = 1 — Dy- The bold dot marks the location where 
the space-time plot in Fig. M has been generated. 



mechanism relies on special properties such as the inter- 
play of several species of particles, competing currents, 
unconventional conservation laws, or special boundary 
conditions. In this paper we have presented physical ar- 
guments why first-order phase transitions are impossi- 
ble in simple two-state models without such attributes 
since short-range interactions between domain walls can- 
not stabilize a fluctuating ordered phase. As an example 
we have revisited the triplet-creation model which was 
believed to display a first-order phase transition. How- 
ever, the first-order behavior is a transient phenomenon 
and eventually crosses over to a continuous phase tran- 
sition, supporting the DP conjecture by Jannsen and 
Grassberger [ p4pq |. Similarly, the annihilation/fission 
process crosses over to a second-order phase transition. 
Thus, first-order phase transitions in one spatial dimen- 
sion are subtle and require a much more robust mecha- 
nism for the elimination of minority islands. 



The temporal evolution of the coarsening process de- 
pends on the dynamics of domain walls between regions 
with low and high particle density. Obviously, high- 
density islands tend to grow until the particle density 
of their environment decreases below a certain thresh- 
old. This threshold, however, depends not only on the 
initial particle density p but also on the total size of all 
active islands. In other words, active islands continue 
to grow until they bind so many particles that the aver- 
age density of particles in their environment is no longer 
sufficient to sustain the growth process. Moreover, large 
islands seem to be more stable than small ones, leading 
to a slow coarsening process. Thus, finite systems evolve 
into a state where only one high-density domain survives. 

The conservation law has another surprising conse- 
quence: Although the system coarsens, the transition is 
still continuous. In fact, plotting the asymptotic station- 
ary value of ps against p — p c , there is no discontinuity 
(see Fig. ||). This can be explained by observing that 
the size of the surviving active domain grows almost lin- 
early with p — p c . Therefore, the reaction-diffusion (|^) 
provides an interesting example of a system with an in- 
stability which exhibits a second-order phase transition 
in one dimension, although the arguments Sect. || cannot 
be applied. 



VI. CONCLUSIONS 



First-order transitions in 1+1-dimensional nonequilib- 
rium models with fluctuating ordered phases require a 
robust mechanism which eliminates islands of the minor- 
ity phase generated by fluctuations. In many cases this 



[2 

[3; 
[4 
[5 

[7; 

[§ 

[9 
[10 

[11 

[12; 
[13 
[14 

[is: 

[16 
[17 



[is; 

[19 



J. Marro and R. Dickman, Nonequilibnum phase tran- 
sitions in lattice models, Cambridge University Press, 
Cambridge, 1999. 

A. L. Toom, in Multicomponent Random Systems, ed. 

by R. L. Dobrushin, in Advances in Probability, Vol. 6 

(Dekker, New York, 1980), pp. 549-575. 

C. H. Bennett and G. Grinstein, Phys. Rev. Lett. 55, 657 

(1985). 

H. Hinrichsen, R. Livi, D. Mukamel, and A. Politi, Phys. 
Rev. E 61, R1032 (2000). 

R. M. Ziff, E. Gulari and Y. Barshad, Phys. Rev. Lett. 
56, 2553 (1986). 

M. R. Evans, Y. Kafri, H. M. Koduvely, and D. Mukamel, 
Phys. Rev. E 58, 2764 (1998). 

C. Godreche, J. M. Luck, M. R. Evans, D. Mukamel, S. 
Sandow, and E. R. Speer, J. Phys. A 28, 6039 (1995). 
K. Oerding, F. van Wijland, J. - P. Leroy, and H. J. Hil- 



horst, eprint cond-mat/9910351 



N. Menyhard and G. Odor, J. Phys. A 28, 4505 (1995). 
N. Menyhard and G. Odor, J. Phys. A 31, 6771 (1998). 
R. Bidaux, N. Boccara, and H. Chate, Phys. Rev. A 39, 
3094 (1989). 

I. Jensen, Phys. rev. A 43, 3187 (1991) 

E. V. Albano, Phys. Rev. Lett. 69, 656 (1992). 

G. Odor, N. Boccara, and G. Szabo, Phys. Rev. E 48, 

3168 (1993). 

G. Odor and A. Szolnoki, Phys. Rev. E 53, 2231 (1996). 
R. Dickman and T. Tome, Phys. Rev. A 44, 4833 ( 1991) 



E. Carlon, M . Henkel, and U. Schollwock, eprint cond 



mat/9912347| to appear in Phys. Rev. E. 

M. J. Howard and U. C. Timber, J. Phys. A 30, 7721 

(1997). 

W. Kinzel, Percolation structures and processes, Ann. 
Isr. Phys. Soc. 5, eds. G. Deutscher and R. Zallen and 
J. Adler, Adam Hilger, Bristol (1983). 



8 



[20] H. Hinrichsen, Nonequilibrium Critical Phenomen a and 



Phase Transi tions into Absorbing States, eprint cond- 



mat/0001070 



[21] P. Grassberger and A. de la Torre, Ann. Phys. (N.Y.) 

122, 373 (1979). 
[22] H. Hinrichsen, A. Jimenez-Dalmaroni, Y. Rozov, and E. 

Domany, Phys. Rev. Lett. 83 4999 (1999); J. Stat. Phys. 

98, 1149 (2000). 
[23] H. Takayasu and A. Y. Tretyakov, Phys. Rev. Lett. 68, 

3060 (1992). 

[24] D. ben Avraham, F. Leyvraz, and S. Redner, Phys. Rev. 

E 50, 1843 (1994). 
[25] J. Cardy and U. C. Tauber, Phys. Rev. Lett. 77, 4870 

(1996). 

[26] N. Menyhard, J. Phys. A 27, 6139 (1994). 

[27] M. H. Kim and H. Park, Phys. Rev. Lett. 73, 2579 (1994). 

[28] H. Hinrichsen, Phys. Rev. E 55, 219 (199 7). 

[29] H. Hinrichsen, eprint :ond-mat/0001177, to appear in 

Phys. Rev. E. 
[30] G. Odor, eprint |cond-mat/0004150[ 
[31] H. Hinrichsen, eprint sond-mat/000434£ . 
[32] B. Kree, B. Schaub, and B. Schmittmann, Phys. Rev. A 

39 (1989), 2214. 
[33] F. van Wijland, K. Oerding, and H. J. Hilhorst, Physica 

A 251, 179. 
[34] H. K. Janssen, Z. Phys. B 42, 151 (1981). 
[35] P. Grassberger, Z. Phys. B 47, 365 (1982). 







